Transfer and bioaccumulation of pesticides in terrestrial arthropods and food webs: State of knowledge and perspectives for research.

Arthropods represent an entry point for pesticide transfers in terrestrial food webs, and pesticide accumulation in upper chain organisms, such as predators can have cascading consequences on ecosystems. However, the mechanisms driving pesticide transfer and bioaccumulation in food webs remain poorly understood. Here we review the literature on pesticide transfers mediated by terrestrial arthropods in food webs. The transfer of pesticides and their potential for bioaccumulation and biomagnification are related to the chemical properties and toxicokinetic of the substances, the resistance and detoxification abilities of the contaminated organisms, as well as by their effects on organisms' life history traits. We further identify four critical areas in which knowledge gain would improve future predictions of pesticides impacts on terrestrial food webs. First, efforts should be made regarding the effects of co-formulants and pesticides mixtures that are currently understudied. Second, progress in the sensitivity of analytical methods would allow the detection of low concentrations of pesticides in small individual arthropods. Quantifying pesticides in arthropods preys, their predators, and arthropods or vertebrates at higher trophic level would bring crucial insights into the bioaccumulation and biomagnification potential of pesticides in real-world terrestrial food webs. Finally, quantifying the influence of the trophic structure and complexity of communities on the transfer of pesticides could address several important sources of variability in bioaccumulation and biomagnification across species and food webs. This narrative review will inspire future studies aiming to quantify pesticide transfers in terrestrial food webs to better capture their ecological consequences in natural and cultivated landscapes.

Biotic Ligand Modeling for the Effect of Major Cations on the Uptake of Cadmium in Folsomia candida Exposed in a Sand-Solution Medium.

Biotic ligand modeling (BLM) approaches are already applied to predict the bioavailability and possible risk of metals in surface water, but need further development for soils. The present study investigated the effect of major cations (Ca2+, Mg2+, Na+, K+, and H+) on cadmium bioaccumulation in the springtail Folsomia candida. To avoid the complexity of real soils and enable control of elemental speciation in the exposure medium, the animals were exposed to different cadmium concentrations in an inert quartz sand-solution medium. Accumulation of cadmium in the animals was measured after 7 days exposure at different cation concentrations. Among the cations, only Ca2+ significantly affected the uptake of cadmium in the springtails. Mg2+ also had higher effects compared with other selected cations. Using a BLM approach, the uptake of cadmium in the animals predicted by taking into account both Ca2+ and Mg2+ activities correlated well with the measured values (R2 = 0.68). The final estimated conditional binding constants for cadmium (log KCd-BL), Ca (log KCa-BL), and Mg (log KMg-BL) of 1.06, 2.14, and 1.23 L/mol, respectively, were in agreement with previously reported values. The match between predicted and measured uptake data confirms the applicability and usefulness of the BLM for predicting the bioavailability of cadmium to springtails and opens the way for its application in soil. Environ Toxicol Chem 2024;43:1090-1096. © 2024 SETAC.

Rainforest transformation reallocates energy from green to brown food webs.

Terrestrial animal biodiversity is increasingly being lost because of land-use change1,2. However, functional and energetic consequences aboveground and belowground and across trophic levels in megadiverse tropical ecosystems remain largely unknown. To fill this gap, we assessed changes in energy fluxes across 'green' aboveground (canopy arthropods and birds) and 'brown' belowground (soil arthropods and earthworms) animal food webs in tropical rainforests and plantations in Sumatra, Indonesia. Our results showed that most of the energy in rainforests is channelled to the belowground animal food web. Oil palm and rubber plantations had similar or, in the case of rubber agroforest, higher total animal energy fluxes compared to rainforest but the key energetic nodes were distinctly different: in rainforest more than 90% of the total animal energy flux was channelled by arthropods in soil and canopy, whereas in plantations more than 50% of the energy was allocated to annelids (earthworms). Land-use change led to a consistent decline in multitrophic energy flux aboveground, whereas belowground food webs responded with reduced energy flux to higher trophic levels, down to -90%, and with shifts from slow (fungal) to fast (bacterial) energy channels and from faeces production towards consumption of soil organic matter. This coincides with previously reported soil carbon stock depletion3. Here we show that well-documented animal biodiversity declines with tropical land-use change4-6 are associated with vast energetic and functional restructuring in food webs across aboveground and belowground ecosystem compartments.

Insights into sequence characteristics and evolutionary history of DGATs in arthropods.

Triacylglycerol (TAG) is crucial in animal energy storage and membrane biogenesis. The conversion of diacylglycerol (DAG) to triacylglycerol (TAG) is catalyzed by diacylglycerol acyltransferase enzymes (DGATs), which are encoded by genes belonging to two distinct gene families. Although arthropods are known to possess DGATs activities and utilize the glycerol-3-phosphate pathway and MAG pathway for TAG biosynthesis, the sequence characterization and evolutionary history of DGATs in arthropods remains unclear. This study aimed to comparatively evaluate genomic analyses of DGATs in 13 arthropod species and 14 outgroup species. We found that arthropods lack SOAT2 genes within the DGAT1 family, while DGAT2, MOGAT3, AWAT1, and AWAT2 were absent from in DGAT2 family. Gene structure and phylogenetic analyses revealed that DGAT1 and DGAT2 genes come from different gene families. The expression patterns of these genes were further analyzed in crustaceans, demonstrating the importance of DGAT1 in TAG biosynthesis. Additionally, we identified the DGAT1 gene in Swimming crab (P. trituberculatus) undergoes a mutually exclusive alternative splicing event in the molt stages. Our newly determined DGAT inventory data provide a more complete scenario and insights into the evolutionary dynamics and functional diversification of DGATs in arthropods.

Exploration and utilization of novel aldoxime, nitrile, and nitro compounds metabolizing enzymes from plants and arthropods.

Aldoxime (R1R2C=NOH) and nitrile (R-C≡N) are nitrogen-containing compounds that are found in species representing all kingdoms of life. The enzymes discovered from the microbial "aldoxime-nitrile" pathway (aldoxime dehydratase, nitrile hydratase, amidase, and nitrilase) have been thoroughly studied because of their industrial importance. Although plants utilize cytochrome P450 monooxygenases to produce aldoxime and nitrile, many biosynthetic pathways are yet to be studied. Cyanogenic millipedes accumulate various nitrile compounds, such as mandelonitrile. However, no such aldoxime- and nitrile-metabolizing enzymes have been identified in millipedes. Here, I review the exploration of novel enzymes from plants and millipedes with characteristics distinct from those of microbial enzymes, the catalysis of industrially useful reactions, and applications of these enzymes for nitrile compound production.

Advances in understanding insect chitin biosynthesis.

Chitin, a natural polymer of N-acetylglucosamine chains, is a principal component of the apical extracellular matrix in arthropods. Chitin microfibrils serve as structural components of natural biocomposites present in the extracellular matrix of a variety of invertebrates including sponges, molluscs, nematodes, fungi and arthropods. In this review, we summarize the frontier advances of insect chitin synthesis. More specifically, we focus on the chitin synthase (CHS), which catalyzes the key biosynthesis step. CHS is also known as an attractive insecticidal target in that this enzyme is absent in mammals, birds or plants. As no insect chitin synthase structure have been reported so far, we review recent studies on glycosyltransferase domain structures derived from fungi and oomycetes, which are conserved in CHS from all species containing chitin. Auxiliary proteins, which coordinate with CHS in chitin biosynthesis and assembly, are also discussed.

Cyclopropane Hydrocarbons from the Springtail Vertagopus sarekensis─A New Class of Cuticular Lipids from Arthropods.

The epicuticle of insects is usually coated with a complex mixture of hydrocarbons, primarily straight-chain and methyl-branched alkanes and alkenes. We were interested in whether springtails (Collembola), a sister class of the insects, also use such compounds. We focused here on Vertagopus sarekensis, an abundant Isotomidae species in European high alpine regions, exhibiting coordinated group behavior and migration. This coordination, suggesting chemical communication, made the species interesting for our study on epicuticular hydrocarbons in springtails with different degrees of group behavior. We isolated a single hydrocarbon from its surface, which is the major epicuticular lipid. The structure was deduced by NMR analysis and GC/MS including derivatization. Total synthesis confirmed the structure as cis,cis-3,4,13,14-bismethylene-24-methyldotriacontane (4, sarekensane). The GC/MS analyses of some other cyclopropane hydrocarbons also synthesized showed the close similarity of both mass spectra and gas chromatographic retention indices of alkenes and cyclopropanes. Therefore, analyses of cuticular alkenes must be performed with appropriate derivatization to distinguish these two types of cuticular hydrocarbons. Sarekensane (4) is the first nonterpenoid cuticular hydrocarbon from Collembola that is biosynthesized via the fatty acid pathway, as are insect hydrocarbons, and contains unprecedented cyclopropane rings in the chain, not previously reported from arthropods.

PEPCK and glucose metabolism homeostasis in arthropods.

The fat body is responsible for a variety of functions related to energy metabolism in arthropods, by controlling the processes of de novo glucose production (gluconeogenesis) and glycogen metabolism. The rate-limiting factor of gluconeogenesis is the enzyme phosphoenolpyruvate carboxykinase (PEPCK), generally considered to be the first committed step in this pathway. Although the study of PEPCK and gluconeogenesis has been for decades restricted to mammalian models, especially focusing on muscle and liver tissue, current research has demonstrated particularities about the regulation of this enzyme in arthropods, and described new functions. This review will focus on arthropod PEPCK, discuss different aspects to PEPCK regulation and function, its general role in the regulation of gluconeogenesis and other pathways. The text also presents our views on potentially important new directions for research involving this enzyme in a variety of metabolic adaptations (e.g. diapause), discussing enzyme isoforms, roles during arthropod embryogenesis, as well as involvement in vector-pathogen interactions, contributing to a better understanding of insect vectors of diseases and their control.

Temporal regulation of temperature tolerances and gene expression in an arctic insect.

Terrestrial arthropods in the Arctic are exposed to highly variable temperatures that frequently reach cold and warm extremes. Yet, ecophysiological studies on arctic insects typically focus on the ability of species to tolerate low temperatures, whereas studies investigating physiological adaptations of species to periodically warm and variable temperatures are few. In this study, we investigated temporal changes in thermal tolerances and the transcriptome in the Greenlandic seed bug Nysius groenlandicus, collected in the field across different times and temperatures in Southern Greenland. We found that plastic changes in heat and cold tolerances occurred rapidly (within hours) and at a daily scale in the field, and that these changes are correlated with diurnal temperature variation. Using RNA sequencing, we provide molecular underpinnings of the rapid adjustments in thermal tolerance across ambient field temperatures and in the laboratory. We show that transcriptional responses are sensitive to daily temperature changes, and days characterized by high temperature variation induced markedly different expression patterns than thermally stable days. Further, genes associated with laboratory-induced heat responses, including expression of heat shock proteins and vitellogenins, were shared across laboratory and field experiments, but induced at time points associated with lower temperatures in the field. Cold stress responses were not manifested at the transcriptomic level.

Polyproline type II helical antifreeze proteins are widespread in Collembola and likely originated over 400 million years ago in the Ordovician Period.

Antifreeze proteins (AFPs) bind to ice crystals to prevent organisms from freezing. A diversity of AFP folds has been found in fish and insects, including alpha helices, globular proteins, and several different beta solenoids. But the variety of AFPs in flightless arthropods, like Collembola, has not yet been adequately assessed. Here, antifreeze activity was shown to be present in 18 of the 22 species of Collembola from cold or temperate zones. Several methods were used to characterize these AFPs, including isolation by ice affinity purification, MALDI mass spectrometry, amino acid composition analysis, tandem mass spectrometry sequencing, transcriptome sequencing, and bioinformatic investigations of sequence databases. All of these AFPs had a high glycine content and were predicted to have the same polyproline type II helical bundle fold, a fold unique to Collembola. These Hexapods arose in the Ordovician Period with the two orders known to produce AFPs diverging around 400 million years ago during the Andean-Saharan Ice Age. Therefore, it is likely that the AFP arose then and persisted in many lineages through the following two ice ages and intervening warm periods, unlike the AFPs of fish which arose independently during the Cenozoic Ice Age beginning ~ 30 million years ago.

CryoEM structure and Alphafold molecular modelling of a novel molluscan hemocyanin.

Hemocyanins are multimeric oxygen transport proteins present in the blood of arthropods and molluscs, containing up to 8 oxygen-binding functional units per monomer. In molluscs, hemocyanins are assembled in decamer 'building blocks' formed of 5 dimer 'plates', routinely forming didecamer or higher-order assemblies with d5 or c5 symmetry. Here we describe the cryoEM structures of the didecamer (20-mer) and tridecamer (30-mer) forms of a novel hemocyanin from the slipper limpet Crepidula fornicata (SLH) at 7.0 and 4.7 Å resolution respectively. We show that two decamers assemble in a 'tail-tail' configuration, forming a partially capped cylinder, with an additional decamer adding on in 'head-tail' configuration to make the tridecamer. Analysis of SLH samples shows substantial heterogeneity, suggesting the presence of many higher-order multimers including tetra- and pentadecamers, formed by successive addition of decamers in head-tail configuration. Retrieval of sequence data for a full-length isoform of SLH enabled the use of Alphafold to produce a molecular model of SLH, which indicated the formation of dimer slabs with high similarity to those found in keyhole limpet hemocyanin. The fit of the molecular model to the cryoEM density was excellent, showing an overall structure where the final two functional units of the subunit (FU-g and FU-h) form the partial cap at one end of the decamer, and permitting analysis of the subunit interfaces governing the assembly of tail-tail and head-tail decamer interactions as well as potential sites for N-glycosylation. Our work contributes to the understanding of higher-order oligomer formation in molluscan hemocyanins and demonstrates the utility of Alphafold for building accurate structural models of large oligomeric proteins.

Tritrophic Interactions among Arthropod Natural Enemies, Herbivores and Plants Considering Volatile Blends at Different Scale Levels.

Herbivore-induced plant volatiles (HIPVs) are released by plants upon damaged or disturbance by phytophagous insects. Plants emit HIPV signals not merely in reaction to tissue damage, but also in response to herbivore salivary secretions, oviposition, and excrement. Although certain volatile chemicals are retained in plant tissues and released rapidly upon damaged, others are synthesized de novo in response to herbivore feeding and emitted not only from damaged tissue but also from nearby by undamaged leaves. HIPVs can be used by predators and parasitoids to locate herbivores at different spatial scales. The HIPV-emitting spatial pattern is dynamic and heterogeneous in nature and influenced by the concentration, chemical makeup, breakdown of the emitted mixes and environmental elements (e.g., turbulence, wind and vegetation) which affect the foraging of biocontrol agents. In addition, sensory capability to detect volatiles and the physical ability to move towards the source were also different between natural enemy individuals. The impacts of HIPVs on arthropod natural enemies have been partially studied at spatial scales, that is why the functions of HIPVs is still subject under much debate. In this review, we summarized the current knowledge and loopholes regarding the role of HIPVs in tritrophic interactions at multiple scale levels. Therefore, we contend that closing these loopholes will make it much easier to use HIPVs for sustainable pest management in agriculture.

Expression of netrin and its receptors uncoordinated-5 and frazzled in arthropods and onychophorans suggests conserved and diverged functions in neuronal pathfinding and synaptogenesis.

BACKGROUND: Development of the nervous system and the correct connection of nerve cells require coordinated axonal pathfinding through an extracellular matrix. Outgrowing axons exhibit directional growth toward or away from external guidance cues such as Netrin. Guidance cues can be detected by growth cones that are located at the end of growing axons through membrane-bound receptors such as Uncoordianted-5 and Frazzled. Binding of Netrin causes reformation of the cytoskeleton and growth of the axon toward (or away from) the source of Netrin production. RESULTS: Here, we investigate the embryonic mRNA expression patterns of netrin genes and their potential receptors, uncoordinated-5 and frazzled in arthropod species that cover all main branches of Arthropoda, that is, Pancrustacea, Myriapoda, and Chelicerata. We also studied the expression patterns in a closely related outgroup species, the onychophoran Euperipatoides kanangrensis, and provide data on expression profiles of these genes in larval tissues of the fly Drosophila melanogaster including the brain and the imaginal disks. CONCLUSION: Our data reveal conserved and diverged aspects of neuronal guidance in Drosophila with respect to the other investigated species and suggest a conserved function in nervous system patterning of the developing appendages.

Characterization of cuticular proteins in CPR family in the wolf spider, Pardosa pseudoannulata, and the response of one subfamily genes to environmental stresses.

Cuticular protein (CP) plays an essential role in the construction and function of exoskeleton in arthropods. CPR family, CP with Rebers and Riddiford (R&R) Consensus, is the largest CP family in insects, but it lacks systematic research in non-insect arthropods. In this study, we explored CPRs in the wolf spider, Pardosa pseudoannulata, a predator to many insect pests. We totally identified 152 CPRs in P. pseudoannulata genome, which were divided into two subgroups based on R&R Consensus sequences, with 12 CPRs in RR-1 and 140 in RR-2. All RR-2 members presented a novel Consensus with 34 amino acids, G-x(8)-G-x(6)-Y-x-A-x(3)-G-x(7)-N-E-x-G, which was a common characteristic for RR-2 CPRs in chelicerates. Transcriptome data was used to document the expression patterns of CPR genes in different tissues and ecdysis processes. The specific expressions were found for part CPR genes, such as five RR-2 genes that were specifically expressed in male genital bulbs and eleven RR-1 genes that were highly expressed in the integument. Due to the limited number and integument-specific expression of RR-1 genes, we further analyzed their responses to different environmental stresses at the transcriptional level. Except for PapsCPR11, ten RR-1 genes responded to at least one environmental stress, among with the expression of PapsCPR12 was significantly changed by three stresses (dryness, low temperature and imidacloprid treatments). Silencing PapsCPR12 increased the tolerance of P. pseudoannulata to imidacloprid. Overall, the results presented novel Consensus characteristics of CPRs in P. pseudoannulata, which was helpful for the identification and evolution analysis of CPRs in non-insect arthropods.

Mitogenome selection in the evolution of key ecological strategies in the ancient hexapod class Collembola.

A longstanding question in evolutionary biology is how natural selection and environmental pressures shape the mitochondrial genomic architectures of organisms. Mitochondria play a pivotal role in cellular respiration and aerobic metabolism, making their genomes functionally highly constrained. Evaluating selective pressures on mitochondrial genes can provide functional and ecological insights into the evolution of organisms. Collembola (springtails) are an ancient hexapod group that includes the oldest terrestrial arthropods in the fossil record, and that are closely associated with soil environments. Of interest is the diversity of habitat stratification preferences (life forms) exhibited by different species within the group. To understand whether signals of positive selection are linked to the evolution of life forms, we analysed 32 published Collembola mitogenomes in a phylomitogenomic framework. We found no evidence that signatures of selection are correlated with the evolution of novel life forms, but rather that mutations have accumulated as a function of time. Our results highlight the importance of nuclear-mitochondrial interactions in the evolution of collembolan life forms and that mitochondrial genomic data should be interpreted with caution, as complex selection signals may complicate evolutionary inferences.

A comprehensive study of arthropod and onychophoran Fox gene expression patterns.

Fox genes represent an evolutionary old class of transcription factor encoding genes that evolved in the last common ancestor of fungi and animals. They represent key-components of multiple gene regulatory networks (GRNs) that are essential for embryonic development. Most of our knowledge about the function of Fox genes comes from vertebrate research, and for arthropods the only comprehensive gene expression analysis is that of the fly Drosophila melanogaster. For other arthropods, only selected Fox genes have been investigated. In this study, we provide the first comprehensive gene expression analysis of arthropod Fox genes including representative species of all main groups of arthropods, Pancrustacea, Myriapoda and Chelicerata. We also provide the first comprehensive analysis of Fox gene expression in an onychophoran species. Our data show that many of the Fox genes likely retained their function during panarthropod evolution highlighting their importance in development. Comparison with published data from other groups of animals shows that this high degree of evolutionary conservation often dates back beyond the last common ancestor of Panarthropoda.

Metabarcoding versus mapping unassembled shotgun reads for identification of prey consumed by arthropod epigeal predators.

BACKGROUND: A central challenge of DNA gut content analysis is to identify prey in a highly degraded DNA community. In this study, we evaluated prey detection using metabarcoding and a method of mapping unassembled shotgun reads (Lazaro). RESULTS: In a mock prey community, metabarcoding did not detect any prey, probably owing to primer choice and/or preferential predator DNA amplification, while Lazaro detected prey with accuracy 43-71%. Gut content analysis of field-collected arthropod epigeal predators (3 ants, 1 dermapteran, and 1 carabid) from agricultural habitats in Brazil (27 samples, 46-273 individuals per sample) revealed that 64% of the prey species detections by either method were not confirmed by melting curve analysis and 87% of the true prey were detected in common. We hypothesized that Lazaro would detect fewer true- and false-positive and more false-negative prey with greater taxonomic resolution than metabarcoding but found that the methods were similar in sensitivity, specificity, false discovery rate, false omission rate, and accuracy. There was a positive correlation between the relative prey DNA concentration in the samples and the number of prey reads detected by Lazaro, while this was inconsistent for metabarcoding. CONCLUSIONS: Metabarcoding and Lazaro had similar, but partially complementary, detection of prey in arthropod predator guts. However, while Lazaro was almost 2× more expensive, the number of reads was related to the amount of prey DNA, suggesting that Lazaro may provide quantitative prey information while metabarcoding did not.

Lack of evidence for conserved parasegmental grooves in arthropods.

In the arthropod model species Drosophila melanogaster, a dipteran fly, segmentation of the anterior-posterior body axis is under control of a hierarchic gene cascade. Segmental boundaries that form morphological grooves are established posteriorly within the segmental expression domain of the segment-polarity gene (SPG) engrailed (en). More important for the development of the fly, however, are the parasegmental boundaries that are established at the interface of en expressing cells and anteriorly adjacent wingless (wg) expressing cells. In Drosophila, both segmental and transient parasegmental grooves form. The latter are positioned anterior to the expression of en. Although the function of the SPGs in establishing and maintaining segmental and parasegmental boundaries is highly conserved among arthropods, parasegmental grooves have only been reported for Drosophila, and a spider (Cupiennius salei). Here, we present new data on en expression, and re-evaluate published data, from four distantly related spiders, including Cupiennius, and a distantly related chelicerate, the harvestman Phalangium opilio. Gene expression analysis of en genes in these animals does not corroborate the presence of parasegmental grooves. Consequently, our data question the general presence of parasegmental grooves in arthropods.

Functionally characterized arthropod pest and pollinator cytochrome P450s associated with xenobiotic metabolism.

The cytochrome P450 family (P450s) of arthropods includes diverse enzymes involved in endogenous essential physiological functions and in the oxidative metabolism of xenobiotics, insecticides and plant allelochemicals. P450s can also establish insecticide selectivity in bees and pollinators. Several arthropod P450s, distributed in different phylogenetic groups, have been associated with xenobiotic metabolism, and some of them have been functionally characterized, using different in vitro and in vivo systems. The purpose of this review is to summarize scientific publications on arthropod P450s from major insect and mite agricultural pests, pollinators and Papilio sp, which have been functionally characterized and shown to metabolize xenobiotics and/or their role (direct or indirect) in pesticide toxicity or resistance has been functionally validated. The phylogenetic relationships among these P450s, the functional systems employed for their characterization and their xenobiotic catalytic properties are presented, in a systematic approach, including critical aspects and limitations. The potential of the primary P450-based metabolic pathway of target and non-target organisms for the development of highly selective insecticides and resistance-breaking formulations may help to improve the efficiency and sustainability of pest control.

Harnessing the functional diversity of plant cystatins to design inhibitor variants highly active against herbivorous arthropod digestive proteases.

Protein engineering approaches have been proposed to improve the inhibitory properties of plant cystatins against herbivorous arthropod digestive proteases, generally involving the site-directed mutagenesis of functionally relevant amino acids or the selection of improved inhibitor variants by phage display approaches. Here, we propose a novel approach where the function-related structural elements of a cystatin are substituted by the corresponding elements of an alternative cystatin. Inhibitory assays were first performed with 20 representative plant cystatins and model Cys proteases, including arthropod proteases, to appreciate the extent of functional variability among the plant cystatin family. The most, and less, potent of these cystatins were then used as 'donors' of structural elements to create hybrids of tomato cystatin SlCYS8 used as a model 'recipient' inhibitor. In brief, inhibitory activities against Cys proteases strongly differed from one plant cystatin to another, with Ki (papain) values diverging by more than 30-fold and inhibitory rates against arthropod proteases varying by up to 50-fold depending on the enzymes assessed. In line with theoretical assumptions from docking models generated for different Cys protease-cystatin combinations, structural element substitutions had a strong impact on the activity of recipient cystatin SlCYS8, positive or negative depending on the basic inhibitory potency of the donor cystatin. Our data confirm the wide variety of cystatin inhibitory profiles among plant taxa. They also demonstrate the usefulness of these proteins as a pool of discrete structural elements for the design of cystatin variants with improved potency against herbivorous pest digestive Cys proteases.